function [x,ix,si]=ps1(A,xmin,xmax,ixmin,ixmax,ov);

if (~exist('ov','var'))
  ov=9;
end

[n,m]=size(A);
assert(n==m);

x=xmin:(xmax-xmin)/ov:xmax;
ix=ixmin:(ixmax-ixmin)/ov:ixmax;

si=zeros(length(ix),length(x));

for k=1:length(x);
  for l=1:length(ix);
    d=svd(A-eye(n)*(x(k)+i*ix(l)));
    si(l,k)=-log(d(n));
    % si(l,k)=svds(A-eye(n)*(x(k)+i*ix(l)),1,0);
  end
end


